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FREE-FALL  DECAY  OF  A  MERCURY  VAPOR  PLASMA  IN  THE  EARLY  AFTERGLOW* 

B.  C.  Gregory 

Physic*  Deportment,  Trent  University,  Peterborough,  Canada 
^  (Received  6  January  1969) 

A  free-fall  theory  la  developed  for  the  decay  of  the  average  electron  number  density 
and  of  the  electron  temperature  In  a  pulsed  cylindrical  mercury  vapor  discharge  plasma 
using  moment  equations  for  the  ions  with  a  collision  term  for  momentum  loss,  and  as¬ 
suming  Maxwellian  electrons.  The  plasma  assumption  wss  used.  Good  agreement  Is  ob¬ 
tained  with  experimental  results  on  four  tubes  of  varying  diameter.1  J 


As  the  mercury  vapor  pressure  is  lowered,  the 
afterglow  decay  profile  of  the  average  electron 
density  < ne(t )>  in  a  cylindrical  discharge  tube  be¬ 
comes  quite  nonexponential  in  time  (Fig.  1 ).  This 
Letter  is  concerned  with  the  nonexponential  be¬ 
havior  of  (,ne(t))  which  occurs  when  the  mean  free 
path  of  the  charged  particles  is  greater  than  the 
transverse  dimensions  of  the  tube.  Experimental 
results  are  explained  by  using  the  first  two  mo¬ 
ment  equations  for  the  ions,  the  plasma  assump¬ 


tion,  an  ion -neutral  momentum -loss  collision 
term,  the  assumption  that  electrons  are  Maxwel¬ 
lian,  and  an  equation  describing  the  energy  decay 
of  the  electrons. 

The  theory  is  based  upon  the  good  agreement 
between  the  simplified  moment  treatment  of  Kino 
and  Shaw*  and  the  more  exact  treatment  of  Self* 
and  Parker.*  We  have  used  the  steady-state  ini¬ 
tial  conditions  for  electron  density  (=  ion  density) 
and  drift  velocity  worked  out  in  Ref.  1  for  cylln- 
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1IG.  1  Measured  average  electron  number-density 
profiles  In  the  afterglow.  Tube  C  with  pressure  as 
parameter . 

drical  geometry,  and  (i)  constant  ionization  with 
no  ion  kinetic  pressure  term,  (ii)  ionization  pro¬ 
portional  to  electron  density  with  no  pressure 
term,  and  (iii)  ionization  proportional  to  electron 
density  with  an  ion  kinetic  pressure  term.  Case 
(iii)  was  shown  in  Ref.  1  to  give  results  very  sim¬ 
ilar  to  the  more  exact  treatment  of  Refs.  2  and  3, 
even  though  the  fluid  method  assumes  that  all  the 
ions  have  the  same  velocity  at  any  point. 

Our  theoretical  results  for  the  afterglow  indi¬ 
cate  that  in  case  (i)  above  the  electron  density  de¬ 
cays  about  20-30%  more  slowly  than  cases  (ii) 
and  (iii),  which  agree  quite  closely.  Since  case 
(iii)  seems  to  be  the  most  realistic  here  we  use 
its  results  for  the  initial  conditions  in  what  fol¬ 
lows. 

In  the  afterglow  we  assume  that  the  ion  and 
electron  densities  are  equal  at  all  points  in  space 
and  at  all  times.  The  ion  kinetic  pressure  tensor 
term  is  dropped.  We  also  assume  that  no  further 
ionization  occurs  after  t  =0.  The  electrons  are 
assumed  Maxwellian: 

ti  -n-n  exp {etpfkT  ),  (1) 

e  eO  e 


where  ne  is  the  electron  number  density,  neQ  lt* 
axial  value,  <p  the  space  potential,  -e  the  elec¬ 
tronic  charge,  k  Boltzmann’s  constant,  and  Te 
the  electron  temperature  (assumed  a  function  of 
time  only  and  not  of  space).  With  these  assump¬ 
tions,  the  firBt  two  moment  equations  for  the  ions 
are 
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and 
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Here  v  is  the  ion  drift  velocity,  v  an  ion-neutral 
collision  frequency,  and  M  the  ion  mass. 

A  third  equation  is  necessary  to  describe  ener¬ 
gy  changes.  This  states  that  the  rate  of  decrease 
of  average  electron  energy  equals  the  electron 
wall  flux  times  the  energy  equivalent  of  the  wall 
potential: 

4:  (a\nkT  2nrdr  =  -2va(vn)  etp  ,  (4) 

dtJo  e  a  w 


where  a  is  the  tube  radius  and  <pw  is  the  tube 
wall  potential.  This  equation  implies  that  each 
electron  arriving  at  the  wall  transports  energy 
etpw  from  the  electron  gas,  which  is  re-Maxwel- 
lianized  immediately.  The  electron  arrival  rate 
at  the  wall  equals  that  of  the  ions.  In  reality 
most  electrons  (energy  less  than  e<pw)  will  under¬ 
go  numerous  wall  collisions  dur'ng  their  life¬ 
times. 

Since  Tg  is  only  a  function  of  time,  Elqs.  (4)  and 
(2)  may  be  combined  to  give 


dT  T  dN 


(5) 


where  N=  J"2irnrdr  and  (3=etp)r/kTe. 

The  fact  that  0  is  almost  constant  for  a  very 
wide  variety  of  discharge  parameters2  enables 
us  to  integrate  numerically  the  three  coupled  non¬ 
linear  equations  (2),  (3),  and  (5)  starting  with  the 
initial  n  and  v  spatial  distributions,  case  (iii) 
above.  Constant  collision  frequency  gives  results 
which  do  not  agree  with  experiment,  and  so  a 
constant  cross  section  for  ion-atom  collisions, 
e,  was  assumed.  The  collision  frequency  is  then 
v  =  KflOf,  where  is  the  atom  density.  A  typical 
set  of  curves  with  kTeQ  =  3,0  eV  and  parameter 
nao  is  shown  in  Fig.  2. 

The  hot-cathode  discharge  tubes  were  described 
earlier.4  Current  pulses  were  of  0.2  to  2.0  A  and 
from  5  to  50  psec  duration.  Most  results  were 
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taken  at  40  Msec,  at  which  time  steady-state  con¬ 
ditions  prevailsd.  The  tubes  were  run  in  an  oven 
at  (100 ±  5)°C  with  a  finger  of  liquid  mercury  pro¬ 
truding  below  into  a  water  bath  whose  tempera¬ 
ture  was  controlled  to  within  ±0.2  C°. 

The  average  electron  number  density  was  mea¬ 
sured  with  a  dipole  resonance  cavity4  and  a  per¬ 
turbation  TMou  mode  cavity.  The  two  methods 
agree  very  well  (see  Fig.  1)  when  care  is  taken 
to  account  for  the  cavity  end  holes  and  the  glass 
tube  envelope.4  The  former  method  was  used 
most  of  the  time  because  the  Q  of  the  dipole - 
mode  cavity  is  higher,  yielding  superior  time 
resolution. 

Electron  temperature  decay  curves  in  these 
rapidly  changing  plasmas  are  exceedingly  diffi¬ 
cult  to  measure.  We  have  used  a  pulsed  Lang¬ 
muir  probe  technique,  a  swept  probe  technique, 
and  a  time- resolving  (1-2  Msec)  3.0-GHz  radiom¬ 
eter.4  However,  at  the  pressures  of  interest 
here  the  absorptivity  of  the  plasma  is  so  low  that 
the  last  method  is  rather  insensitive.  We  chose 
to  use  the  radiometer  to  verify  the  probe  results 
at  higher  pressures.  Even  so,  we  do  not  place 
much  confidence  in  our  probe  results  at  very  low 
pressures. 


FIG.  2.  Theoretical  normalized  average  electron 
number  density  Parameter  njf.  *T»o*30eV. 


Electron  density  profiles  were  extremely  re¬ 
producible,  even  after  periods  of  months.  The 
shape  of  the  profile  did  not  depend  on  the  dis¬ 
charge  current  pulse,  but  only  on  pressure  and 
initial  electron  temperature.  For  tubes  of  differ¬ 
ent  radii  but  with  roughly  similar  initial  condi¬ 
tions,  approximate  time  constants  of  the  (ne(t)) 
profiles  scaled  directly  as  the  tube  radii,  as 
might  bo  expected  for  a  free-fall  regime. 

To  compare  theory  and  experiment  we  proceed 
as  follows:  (i)  Choose  an  initial  electron  temper¬ 
ature  (usually  measured)  and  select  initial  densi¬ 
ty  and  velocity  space  profiles.  We  have  used 
those  corresponding  to  free-fall  conditions  as  de¬ 
scribed  above,  but  collision  effects  could  be  ac¬ 
counted  for  here  using  the  theories  of  Self  and 
Ewald7  and  Forrest  and  Franklin.4  (ii)  Calculate 
naa.  Alternatively  this  may  be  selected  to  best 
fit  experimental  data,  (iii)  Solve  Eqs.  (2),  (3), 
and  (5)  to  obtain  (»(/)>  and  kTe(t). 

Figure  3  shows  such  a  comparison  for  two 
tubes.  Agreement  for  (n(f)>  is  excellent,  but  for 
kTe(t)  somewhat  poorer.  As  stated  above,  we 
are  not  very  confident  of  the  experimental  elec¬ 
tron  temperature  points,  especially  later  in  the 
afterglow. 


FIG.  3.  Comparison  of  experimental  electron  densi¬ 
ty  and  temperature  profiles  in  the  afterglow  for  two 
tubes  (radii  1.13  and  1.85  cm). 
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The  parameter  on^  used  to  best  fit  theory  with 
experiment  yields  <r*  5*  10“ 10  m*.  a  value  of  the 
order  of  the  elastic  ion-atom  cross  section  but 
at  least  an  order  of  magnitude  smaller  than  the 
average  charge-transfer  cross  section.*  This  is 
not  yet  understood  fully.  Possibly  the  introduc¬ 
tion  of  the  ion  kinetic  pressure  term  in  the  after¬ 
glow  theory  will  clarify  this  point. 

Using  this  theory  with  experimentally  deter¬ 
mined  initial  electron  temperatures,  we  were 
able  to  obtain  excellent  agreement  with  electron 
density  measurements  made  on  four  discharge 
tubes  of  radii  0.55,  0.85,  1.13,  and  1.85  cm. 
More  complete  results  will  be  presented  else¬ 
where.  This  work  tends  to  support  the  validity 
of  using  fluid  theory  for  this  sort  of  problem.  It 
also  predicts  supercooling  of  the  electrons10; 
whether  or  not  this  is  real  or  a  phenomenon  in¬ 
troduced  by  approximations  in  the  theory  will  be 
explored. 
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